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Abstract. After a brief review of the Minimal Supersymmetric Standard Model 
(MSSM) and specifically the Minimal Supergravity Model (SUGRA), the prospects 
for discovering and studying SUSY at the CERN Large Hadron Collider are reviewed. 
The possible role for a future Lepton Collider — whether [i + or e+e~ - is also 
discussed. 



I INTRODUCTION 

The many attractive features of the Minimal Supersymmetric Standard Model [1] 
or MSSM have made it a leading candidate for physics beyond the Standard Model. 
Of course there is no direct evidence for SUSY. The current limits [2] on SUSY 
masses from LEP are close to its ultimate kinematic reach. LEP will extend the 
limits on a Higgs boson from the present 77GeV [2] up to ^ 95 GeV [3]. Discovery 
of a light Higgs would not prove the existence of SUSY but would be a strong 
hint: the light Higgs boson must have a mass less than 130 GeV in the MSSM 
and less than 150 GeV in a rather general class of SUSY models [4], while in the 
Standard Model it must be heavier than about 130 GeV if the theory holds up to 
a high scale [5]. The next run of the Tevatron will have a better chance to find 
SUSY particles; the channel X1X2 ~^ ^ + ^~Xi^ ±iy Xi can De sensitive to masses up 
to ~ 200 GeV for some choices of the other parameters [6]. But the reach in this 
channel is quite model dependent. 

The definitive search for weak-scale SUSY, therefore, will have to await the LHC. 
The LHC with 10 fb" 1 , 10% of its design luminosity per year, can detect gluinos 
and squarks up to about 2 TeV in the multi-jet plus missing transverse energy $t 
channels [7-9] compared to an expected mass scale of less than 1 TeV [10]. It 
is difficult to reconstruct masses directly because every SUSY event contains two 
missing lightest SUSY particles Xv It is possible, however, to use endpoints of 
kinematic distributions to determine combinations of masses. In favorable cases 
these combinations can be used in a global fit to determine the model parameters. 
If SUSY is indeed the right answer, there should be a lot known about it before 



the Next Lepton Collider — whether fi + fi~ or e + e~ - is built. It is probably 
difficult to study the whole SUSY spectrum at the LHC, however, so an NLC is 
also expected to play an important role. 

II MINIMAL SUSY STANDARD MODEL 

The Minimal Supersymmetric Standard Model [1] (MSSM) has for each Standard 
Model particle a partner differing in spin by A J = 1/2. For each gauge boson there 
is a J = 1/2 gaugino, and for each chiral fermion there is a scalar sfermion. Two 
Higgs doublets and their corresponding Higgsinos are needed to give masses to all 
the quarks and to cancel anomalies. The SUSY particles have couplings determined 
by supersymmetry and are degenerate in mass with their Standard Model partners. 

There is at present no experimental evidence for SUSY. There is one possible 
experimental hint: the renormalization group equations imply that the SU(3) x 
SU(2) x U(l) gauge couplings measured at the Z mass meet in a way consistent 
with grand unification for the MSSM with Msusy ~ 1 TeV but not for the Standard 
Model. [12] The unification is actually not quite perfect, but it is within the range 
that could be covered by GUT threshold corrections. Other possible hints have 
been widely discussed but have mostly been discredited. 

SUSY must of course be broken, since there is certainly no selectron degen- 
erate with the electron. It is not possible to obtain an acceptable spectrum by 
breaking SUSY spontaneously using only the MSSM fields. However, mass terms 
for gauginos, Higgsinos, and sfermions do not break the £77(3) x SU(2) x U(l) 
gauge invariance of the MSSM, so they can be added by hand without spoiling its 
renormalizability. There are also soft bilinear (B) and trilinear (Aijk) couplings 
that are also gauge invariant and can be added. Finally, it is necessary to add a 
SUSY-conserving Higgsino mass (//). It is generally assumed that SUSY breaking 
occurs spontaneously in a "hidden sector" and is communicated to the MSSM via 
some common interaction such as gravity. If SUSY is discovered, understanding 
the mechanism for its breaking will become an important issue in particle physics. 

After SUSY is broken, all the states with the same quantum numbers mix. The 
7, Z, Hi, andH 2 mix to give four neutralinos X? 2 3 4- The W ± andH ± mix to give 
two charginos xf 2 - The left and right squarks and sleptons also mix; this mixing is 
proportional to the fermion mass and so is significant only for the third generation. 

The most general MSSM allows baryon and lepton number violation, giving 
proton decay at the weak scale. The simplest solution is to impose invariance 
under a discrete symmetry 

R = (_1)3(B-L)+2S 

Note that R = +1 for all Standard Model particles, and R = -1 for all SUSY 
particles. Thus R parity invariance implies that SUSY particles are produced in 
pairs, that they decay to other SUSY particles, and that the lightest SUSY particle 
(LSP) is absolutely stable. Cosmological constraints then require that the LSP be 
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FIGURE 1. Evolution of masses in the SUGRA model, from Ref. 16. 

neutral and weakly interacting, so that it escapes from any detector. Thus, the 
basic signature of (R parity conserving) SUSY is missing transverse energy fx- 

Ill MINIMAL SUGRA MODEL 

The most general MSSM has more than 100 parameters. Many recent phe- 
nomenological studies have been based on a more restrictive model, the mini- 
mal supergravity (SUGRA) model [13]. This is very similar to the Constrained 
MSSM [14], although the latter adds some additional constraints. If SUSY break- 
ing is communicated through gravity, then it is plausible that the SUSY breaking 
terms, like gravity, are universal at the GUT scale. In particular, if all the scalar 
masses are identical, then electroweak symmetry must be unbroken. It turns out 
that the Clebsch-Gordon coefficients are such that when the parameters are run 
from the GUT scale to the weak scale using the renormalization group equations, 
the large top Yukawa coupling drives the mass-squared of the Higgs negative, break- 
ing electroweak symmetry but not color or charge. A example of this evolution is 
shown in Fig. 1. It is convenient to eliminate B and /z 2 in favor of Mz and tan/?. 
Then the minimal SUGRA model is characterized by just four parameters and the 
sign of \x: 

• mo: the common squark, slepton, and Higgs mass at Mqut- 

• rax j2- the common gaugino mass at Mqut- 

• A : the common trilinear coupling at Mqut- 

• tan/? = v\jvi'- the ratio of Higgs vacuum expectation values at Mz- 




• sgn fj, — ±1. 




500 1000 1500 2000 500 1000 1500 2000 
m (GeV) m (GeV) 

FIGURE 2. SUGRA discovery limits in the SUGRA model at the LHC with lOftT 1 integrated 

luminosity. 0£: ft? + jets + no leptons. 11: fix + jets + one lepton. OS: Opposite sign dileptons. 

SS: same sign dileptons. 3£: trileptons. 2£,0j: dileptons with jet veto (from slepton production). 

31, Oj: trileptons with jet veto (from gaugino production). From Ref. 11. 

It turns out that Aq is not very important for phenomenology at the weak scale. 

While the SUGRA model provides a much more tractable parameter space than 
the general MSSM, one should remember that it is only one possible model. It is 
possible that the assumption of universal masses is not correct. It is also possible 
that SUSY breaking is communicated at a much lower mass scale, as in gauge 
mediated models [15]. 

IV SUSY SIGNATURES AT LHC 

The LHC is a pp collider with an energy of 14 TeV, enough to produce g and q 
pairs with ^ 2 TeV even with 10 fb" 1 . These are typically produced with px ~ M, 
so they move slowly in the lab frame and their decay products are widely separated. 
If R parity is conserved, they will decay into the LSP Xi pl us multiple jets and 
perhaps multiple leptons. A typical decay chain might be: 

9 -> Ql + q 



Such decay chains produce many possible signatures combining jets, leptons, and 
]pT from Xi an d ^'s. For example, since the g is self-conjugate, there are isolated 
same-sign dileptons £ ± £ ± . 

The 5<T discovery limits at the LHC for 10 ftT 1 integrated luminosity are shown 
in Figure 2. Note that the reach is > 2 TeV in the missing energy channels and > 
1 TeV in the multi-lepton channels. Thus the LHC should find multiple signatures 
for SUSY with only 10 ftT 1 if it exists at the weak scale. 

V PRECISION MEASUREMENTS AT LHC 

While it is easy to find signals for SUSY at LHC, there are two missing x5' s m 
each event, making it difficult to reconstruct masses. However, it is possible [18] to 
exploit the cascade decays characteristic of the MSSM to determine combinations 
of masses. The strategy is to start at the bottom of the decay chain and work 
up, partially reconstruct specific final states and relating precision measurements 
of endpoints of kinematic distributions to combinations of masses. A global fit to 
these combinations can then be used to determine the model parameters, at least in 
favorable cases. It would be better to make the global fit not just to such endpoints 
but to all distributions, but this is more difficult technically and perhaps premature 
at this stage. 

What combinations of masses can be determined in this way depends on the decay 
modes and so requires study of specific SUSY models. The CERN LHC Program 
Committee (LHCC) chose five SUGRA points for detailed study by the ATLAS 
and CMS Collaborations. The parameters of these points and some representative 
masses are listed in Table 1. Point 3 is the "comparison point," selected so that 
every existing or proposed accelerator can discover something. Point 5 was chosen 
to give the right cold dark matter for cosmology and so is perhaps the most realistic. 
Points 1 and 2 have gluino and squark masses of order 1 TeV. Point 4 has the 
squarks much heavier than the gluinos. It is close to the boundary allowed by 
electroweak symmetry breaking (at least with ISAJET 7.22) and so has a relatively 
small \i and a large mixing between gauginos and Higgsinos. Studying these specific 
points has proved surprisingly useful [17-20]. 

VI EFFECTIVE MASS 

Gluinos and squarks are strongly produced at the LHC. But production cross 
sections fall rapidly with the produced mass, so it is important to find a variable 
that measures the produced mass for events with missing particles. A variable that 
works well for SUSY is the effective mass, defined as the sum of the missing energy 
and the prs of the first four jets, 



M eS = $ T + pr A + p Tj2 + PT,3 + PT,4 



TABLE 1. Parameters of the LHCC SUGRA points and some representative masses from 
ISAJET 7.22 [21]. 
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Samples of signal and Standard Model background events were generated with 
ISAJET [21]. To separate SUSY from the Standard Model background, events 
were selected with multiple jets plus missing energy: 

• $ T > max(100GeV,0.2M cfr ); 

• > 4 jets with p T > 50 GeV and p Tjl > 100 GeV; 

• Transverse sphericity St > 0.2; 

• No fi or isolated e with pt > 20 GeV, 77 < 2.5. 

Then the SUSY signal emerges from the Standard Model background for large M e s, 
as can be seen from Figure 3. The signal cross section is of order 1 pb in the region 
where it dominates, so it could be discovered in about one month at 10 32 cm -2 s -1 . 
(Of course, it would take much longer than this to understand the detectors.) 

The value of M b q at which the signal emerges from the background scales with 
the SUSY mass scale. To test this, 100 random SUGRA models were generated, 
and the peak of the M cS signal was compared with the SUSY mass scale, defined 
by 

M SUSY = mm{M § , M fi ) 

The scatter plot, shown in Fig. 4, shows a good correlation between the peak and 
the SUSY mass scale, allowing one to determine the SUSY mass scale to about 
10%. 

VII RECONSTRUCTION OF SPECIFIC FINAL 

STATES 

The precision measurements of specific combinations of masses are based on the 
partial reconstruction of the corresponding final states. For each case, SUSY and 
Standard Model background event samples were generated with ISAJET [21] or 
PYTHIA [22], a simple particle- level detector simulation incorporating resolutions 
characteristic of ATLAS and CMS was made, and analysis cuts as described below 
were applied. 
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FIGURE 3. M e g distributions after cuts. Open circles: SUSY signal. Solid circles: it. Upward 
triangles: W + jets. Downward triangles: Z + jets. Squares: QCD jets. Shaded histogram: Sum 
of Standard Model backgrounds. From Ref 18. 
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FIGURE 4. Scatter plot of signal peak in M e g vs. Msusy defined in text. From Ref. 18. 



A Measurement of M(%°) - M{%\) 



The prototype of all the precision measurements is based on the decay X2 
Xi£ + £~ at Point 3. Point 3 has unusual branching ratios: 



B(g 616 + h.c.) = 89% 
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FIGURE 5. distribution for SUSY events at Point 3 (histogram) and for Standard Model 

background (shaded) after cuts described in the text. From Ref .18. 

B{bx - x° 2 b) = 86% 
B(x°2 - xTn = 2 x 17% 

The dominant decay of g — * bj) arises because the b\ is lighter than the g but the 
other squarks are heavier. Events were selected to have two leptons and two b jets: 

• £ + £~ pair with p T)t > 10 GeV, i] < 2.5. 

• > 2 jets tagged as b quarks with p T > 15 GeV and rj < 2. 

• No $t cut was used. 

All distributions shown include a 60% tagging efficiency for 6's and 90% efficiency 
for leptons within the kinematic cuts given above. 

The result of this analysis is a spectacular edge at M(x°) — M(xi) endpoint with 
almost no Standard Model background, as can be seen in Figure 5. Most of the 
SUSY background comes from two xt decays and can be removed by plotting the 
distribution for 

e + e~ + — 2e ± /f f 

This analysis clearly would have huge statistics and would be much easier than 
measuring M\y at Tevatron. Given the current M\y results, it seems conservative 
to estimate an error 

A(M(x!D-M(x?))=50MeV 

for 10 fb _1 . 

Point 3 is perhaps unusually easy, but there is a similar edge at Point 4 plus Z 
peak from heavier gauginos, as can be seen from Figure 6. The estimated error in 
this case is ±1 GeV. 
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FIGURE 6. £+£- distribution for SUSY events at Point 4. From Ref .23. 




FIGURE 7. Projections of the gluino-sbottom mass scatter plot. From 18. 



B g and 61 Reconstruction 

The next step at Point 3 is to combine an £ + £~ pair near the edge with jets to 
determine the b\ and g masses. Events were selected with 

• > 2 jets tagged as b jets with pr > 15 GeV, r\ < 2; 

• £ + £- pair with 45 < M(££) < 55 GeV. 



For an £ + £ pair near the endpoint, the Xi must be soft in the x° res t frame, so 
that 

where M(xi) must be determined from a global fit. The approximately recon- 
structed X2 was combined with one of masses coming from combining the x\ with 
one b to make M(b\) and with a second b to make M(g) using the correct Xi mass. 
The resulting projections, shown in Figure 7, display good resolution on the mass 
difference between the g and the b\ masses — just like for D* — > D7r. By varying 
the assumed x? mass, one finds 

AM(6i) = ±1.5AM(x?) ± 3 GeV 
A (m(£) - M (61)) = ±2 GeV 

C Reconstruction of h — > 66 

For Point 5, the decay X2 ~~ * Xi^ is kinematically allowed and has a branching 
ratio of 64%. Events were selected with 

• > 4 jets with p T > 50 GeV, p Tjl > 100 GeV; 

• Transverse sphericity St > 0.2; 

• M cff = $ T + E ti Pr,i > 800 GeV; 

• $ T > max(100GeV,0.2M cff ). 

and the bb mass was plotted for all pairs of b jets with p T b > 25 GeV and r] b < 2. A 
correction factor was applied to the measured b jet energies to account for neutrinos 
and energy loss out of the cone, and a 60% 6-tagging efficiency was assumed. The 
resulting distribution, Figure 8, has a peak at the Higgs mass with a substantial 
SUSY background but very little Standard Model background. This signal is much 
easier than h — > 77 and would be the discovery mode for the Higgs at this point. 
The h — > bb candidates can be used to reconstruct the decay chain 

9 + 9 ^QLq + QrQ 
Ql -> X2Q -> Xihq, qn -> X°iQ 

To select these events exactly two additional jets with pr > 75 GeV were required. 
Then since one of the two qbb combinations comes from the squark decay, the 
smaller of them must have an endpoint at a function of the squark mass and the 
other masses in the problem. The squark mass can be measured to about 40 GeV 
in this way. 
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FIGURE 8. M(bb) for pairs of b jets for the Point 5 signal (open histogram) and for the sum 
of all backgrounds (shaded histogram) . The smooth curve is a Gaussian plus quadratic fit to the 
signal. The light Higgs mass is 104.15 GeV. 



D £+£- Again 



Consider dileptons for Point 5. The mass distribution after the by now standard 
cuts shows a dramatic edge at about 109 GeV. Since this decay must compete with 
the two-body decay Xo. ~^ Xih, it cannot be a direct three-body decay X2 Xi^ + ^~- 
In fact it comes from two sequential two-body decays, X2 ~* — > Vi^ ± ^ =F , and 
the edge determines 



M max (££) = M{x 



M| 
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to about 1 GeV. One should do a complete fit to the Higgs and dilepton events to 
extract the maximum information on all the masses. This has not yet been done. 
The variable most sensitive to the slepton mass is Pt,2/pt,i- This distribution was 
compared for two different values of m , from which it seems that the slepton mass 
can be estimated to AM(4) ~ 3 GeV. 



E Measurement of M(g) - M(%°), M(xf) 

Gluinos dominate at Point 4 since tuq is large. This means that there is a lot of 
combinatorial background from the many jets in the final state. The strategy of 
this analysis [19,23] is to use trilepton events to select the process 



9 + 9 -> xlw + XiQQ 



X2->iiJ->ix° 
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FIGURE 9. Mu for the Point 5 signal (open histogram) and the sum of all backgrounds (shaded 
histogram). From Ref. 24. 
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FIGURE 10. Dijet mass distributions for Point 4 after cuts. The dashed curve shows only the 
right pairing based on generator information. From Ref. 23. 



Then the dijet mass distributions for the right jet pairing have a common endpoint 
since M(*°) « M{xf). 

Events were selected by requiring three leptons and four jets: 

• 3 isolated L p T > 20, 10, 10 GeV, \rj\ < 2.5 

• One opposite-sign, same-flavor lepton pair with Mu < 72 GeV. 

• 4 jets; p T > 150, 120, 70, 40 GeV, \r]\ < 3.2. 
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FIGURE 11. Errors from global fits to the SUGRA parameters. Circles: Fit I. Squares: Ulti- 
mate Fit II. 



• No additional jets with px > 40 GeV and \t]\ < 5 to minimize combinatorics. 

No $t cut was used. With these cuts there are 250 signal events, 30 gq background 
events, and 18 other SUSY and Standard Model background events for 30 fb~ . The 
pairing between the two highest and the two lowest px jets is usually wrong and 
so was eliminated. The dijet mass distribution for the other pairings is shown in 
Figure 10. There is an endpoint at about the right point. 



VIII FITTING SUGRA PARAMETERS 

The precision measurements described here are only a fraction of those in 
Refs. 18, 23, and 24. Ideally one should combine these measurements with a large 
number of other ones and do a global fit, but this would require generating many 
signal signal event samples. Instead a much simpler procedure has been adopted. 
SUGRA parameters were generated at random, the mass spectrum for each SUGRA 
point was calculated, and these were compared with the precision measurements 
and their estimated errors. 

Two such fits have been made. Fit I [18] uses only the measurements developed in 
Ref. 18. It assumes that the Higgs mass can be related to the SUGRA parameters 
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FIGURE 12. Solid: Visible rr mass distribution for 3-prong r decays. Dashed: Contribution 
from X2 i" 7 - Dash-dotted: Contribution from higher mass gauginos. 

with an error AMh = 3 GeV, about the current theoretical error. It bases the 
statistical errors on 10 fb _1 . Fit II [25] adds some additional precision measurements 
developed after Ref. 18. It adds some other data, e.g., from changing the squark 
masses at Points 1 and 2 and seeing the effect of this on the mean of the 
hardest jet. This is not fully justified, but it is a plausible way of estimating 
the improvement from fitting some of the kinematic distributions as well as the 
precision measurements. The ultimate version of Fit II also assumes a theoretical 
error on the Higgs mass less than the expected experimental error from h — > 77, 
AMh = 0.2 GeV and scales the statistical errors to 300fb _1 . 

Both fits scanned the SUGRA parameter space and determined the 68% confi- 
dence interval for each parameter. The results are summarized in Figure 11. No 
disconnected solutions were found. In particular, sgn fi was correctly determined, 
although this required including additional information in the fit in some cases. 
The gluino and squark masses are insensitive to mo at Points 1 and 2, and there 
is no information available on the slepton masses, accounting for the larger errors 
on m at these points. Finally, A is poorly constrained in all cases, even for the 
ultimate version of Fit II. It is possible to determine the weak scale parameters A t 
and Ab, but these are insensitive to A . 

IX r MODES AT LARGE tan/3 

The five LHCC points do not exhaust the possibilities of even the minimal 
SUGRA model. For example, the f\ is light for large tan/5, so r decays can be 
dominant. Consider the SUGRA point m = = 200 GeV, A Q = 0, tan (3 = 45, 



H < 0. Then 

B(X°2 - ^M) » ^(xf - ^V) » 100% 

Discovery of the SUSY signal is still straightforward, but none of the precision 
measurements discussed above are applicable. 

One possible approach is to require two 3-prong hadronic r's to maximize the 
visible rr mass and hence its sensitivity to the endpoint analogous to that dis- 
cussed in Section VII. D. The difference of t + t~ and t ± t ± is used to eliminate the 
contribution from two Xi decays. The resulting distribution is shown in Figure 12. 
There is clearly an endpoint visible from the contribution of X2 ~^ plus a con- 
tribution from heavier gauginos. Signatures like this require more study both for 
the Tevatron and for the LHC. 

X OUTLOOK FOR LEPTON COLLIDERS 

If SUSY exists at the electroweak scale, it should be straightforward to find 
signals for it at the LHC. It is possible in many cases to make precision measure- 
ments, and if the SUSY model is relatively simple, these can be used to determine 
its parameters. 

The LHC will mainly produce gluinos and squarks. In SUGRA these tend to 
decay mainly into the lighter gauginos; the heavier ones are dominantly Higgsino 
and so are suppressed both by their masses and by their couplings. The direct 
production of sleptons and sneutrinos is also very small, although they can be 
produced as decay products of the light gauginos if they are light enough. Finally, 
the heavy Higgs bosons have small production cross sections, and their dominant 
decay modes have large backgrounds. One should not underestimate the ingenuity 
of experimentalists with real data, but it seems likely that the LHC will not be 
able to study the entire SUSY spectrum. 

A Next Lepton Collider with ~ 500 GeV should be able to detect any SUSY 
particles except the Xi that are kinematically accessible. Sleptons probably rep- 
resent the best opportunity to make significant progress beyond what has been 
learned from the LHC. One of the attractive features of the SUGRA model is that 
the Xi is a good dark matter candidate, and the abundance of cold dark matter 
favors light sleptons [26]. A lepton collider provides an important additional con- 
straint that the slepton pairs are produced with known energy, and this allows 
precise measurements to be made [27,28]. But if more than one slepton is being 
produced, the spectrum can be quite complex, so high luminosity (as well as enough 
energy) may be essential. 
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